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Abstract. In comparison with other well studied star formation regions, Taurus is unusual in several respects, (i) Its stellar 
initial mass function (IMF) peaks at relatively high mass (~ 0.8M G ), but contains very few stars much more massive than 1M Q , 
and is relatively deficient in brown dwarfs, (ii) It has a higher binary fraction, particularly at large separations, (iii) Its core mass 
function is strongly peaked at a few M Q , and the cores have extended envelopes and relatively low levels of turbulence. 

We present here the results of an ensemble of hydrodynamic simulations which suggest that the unusual stellar IMF in Taurus 
is a direct consequence of the unusual properties of the cores there. By following the collapse and fragmentation of cores having 
properties typical of Taurus, we find that roughly 50% of the objects formed in a core, predominantly the low-mass ones, are 
ejected from the core to form a population of low-mass stars and brown dwarfs with a flat mass function. The remaining objects 
form multiple systems within the core, accreting until their masses approach 1M S ; this produces a population of intermediate- 
mass stars whose mass function peaks at ~ 0.8M o . Together these two populations reproduce the IMF in Taurus very well. 
This demonstrates, for the first time, a direct causal link between the core mass function and the stellar IMF in a star formation 
region. 

Key words. Stars: formation - ISM: clouds-structure - Methods: numerical 



1. Introduction 

Stars form within dense molecular cores (e.g. Andre et al., 
2000), and the densest and most centrally condensed cores, i.e. 
those closest to forming stars, are known as prestellar cores 
(Ward-Thompson et al., 1994). In Ophiuchus (Motte, Andre & 
Neri, 1998), Serpens (Testi & Sargent, 1998) and Orion (Motte 
et al., 2001), the mass function of pre-stellar cores is remark- 
ably similar to the IMF for field stars and clustered star for- 
mation regions. This suggests that there is a simple mapping 
from the core mass function into the stellar IMF, with the mean 
masses of the stars forming within a core being proportional to 
the core mass. However, the details of how this works are still 
uncertain. 

In this paper we present numerical simulations of core col- 
lapse and fragmentation which demonstrate a causal link be- 
tween the core mass function and the stellar IMF in the Taurus 
star formation region. We choose to study the Taurus region, 
because it is nearby and extended on the sky, and has therefore 
been studied in detail. Molecular-line mapping has yielded es- 
timates of the Taurus core mass function (Onishi et al., 2002), 
and deep surveys of its stellar content have revealed the Taurus 
IMF down to the deuterium burning limit (Luhman et al., 
2003a). 
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2. Observational background 

Determining and explaining the stellar initial mass function 
(IMF) - i.e. the number of stars, Nm dM, forming with mass 
in the interval (M, M + dM) - is one of the longest standing 
problems in the study of star formation. 

Salpeter (1955) first investigated the IMF for field stars in 
the solar neighbourhood and found a power-law fit, Nm x M~ a 
with a ot 2.35, for stars in the mass range (O.4M < M < 
1OM ). Subsequent work (e.g. Miller & Scalo, 1979; Scalo, 
1986; Kroupa, 2002) has extended the range to both lower and 
higher masses. In the field, Kroupa (2002) finds that Salpeter's 
result is still a good approximation for intermediate and high 
masses (M > O.5M ), but at lower masses the stellar IMF flat- 
tens (i.e. a is smaller). Specifically, a ~ 1.3 + 0.5 for low-mass 
stars, 0.08M Q < M < 0.5M Q ; and a ~ 0.3 + 0.7 for brown 
dwarfs, O.O1M < M < O.O8M . 

In individual star formation regions, there appears to be 
some variance in the stellar IMF, but in regions of clustered 
star formation this variance is small. For example, Muench et 
al. (2002) and Luhman et al. (2003b) find stellar IMFs broadly 
similar to Kroupa's field star IMF in - respectively - Orion and 
IC348. 

However, in Taurus, which is a region of distributed star 
formation, the stellar IMF is markedly different. Figurenjcom- 
pares the IMF of Taurus (solid-line histogram, adapted from 
Luhman et al., 2003a) with that of Orion (shaded histogram, 
adapted from Muench et al., 2002). We note the following 
points from Figure [2 
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Fig. 1. A plot of number of stars, N, versus logarithmic stellar mass, togio\M\, comparing the IMF in Taurus (solid-line histogram, 
adapted from Luhman et al., 2003a) with that in Orion (shaded histogram, adapted from Muench et al., 2002). The Orion IMF is 
normalised to have the same number of stars as the Taurus IMF. The two IMFs are clearly different. 



- The Taurus IMF peaks at a higher mass, ~ O.8M (as com- 
pared with ~ O.2M in Orion). 

- The Taurus IMF contains very few stars more massive than 
1M (whereas Orion contains several stars above 10M o ). 

- The Taurus IMF contains many fewer brown dwarfs per star 
than Orion (less than half as many) 

- At the low-mass end the Taurus IMF is approximately flat. 

There is further evidence that the pattern of star formation 
in Taurus is different from that in regions of clustered star for- 
mation. 

- There is a higher binary fraction in Taurus than in Orion, 
particularly at large separations. 

- The prestellar cores in Taurus (Onishi et al., 2002) have 
a narrower mass function than those in Orion (Motte et 
al., 2001) or Ophiuchus (Motte et al., 1998); most cores 
in Taurus have M ~ 5 ± 3M Q . 

- The prestellar cores in Taurus are on average more extended 
and less turbulent than those in Orion and Ophiuchus. 

3. Modelling the collapse and fragmentation of 
pre-stellar cores in Taurus 

3.1. Initial conditions and constitutive physics 

In view of the strongly peaked core mass function in Taurus, 
we consider a single core mass, M = 5M . The initial condi- 
tions are chosen to mimic the cores observed in Taurus (Ward- 
Thompson et al., 1994; Ward-Thompson, Motte & Andre, 
1999; Jijina, Myers & Adams, 1999; Ward-Thompson, Andre 
& Kirk, 2002; Kirk, Ward-Thompson & Andre, 2003). The 
core density profile is initially 
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with pkemei = 3 x 10 18 gcm 3 (corresponding to n H , - 7 X 
10 5 cirr 3 ), r kerae i = 5, 000 AU and r boundary = 50, 000 AU. 

The isothermal sound speed in a core is initially co = 
0.19 km s~' (corresponding to molecular gas at 10 K), and so 
the initial Jeans mass in the centre of a core is ~ O.8M , and 
the initial ratio of thermal to gravitational potential energy is 
^thermal/ 1^1 = 0.45. In order to capture the switch from isother- 
mality to adiabaticity which is expected to occur as the density 
rises above p crrt ~ 10~ 13 gcrrT 3 (eg. Tohline, 1982; Masunaga 
& Inutsuka, 2000), we use the equation of state 
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A divergence-free Gaussian random velocity field is super- 
imposed on each core to simulate turbulence. The power spec- 
trum of the turbulence is P(k) oc kr A . Fifty realizations of this 
representative core are treated. The only difference between the 
different realizations is the random number seed used to gener- 
ate the turbulent velocity field, and the initial ratio of turbulent 
to gravitational potential energy, t/turbuient/PI- ln order to match 
the levels of turbulence measured in Taurus, twenty cores have 
£/turbuient/|£2| = 0.05, twenty cores have f/ tur buient/l^l = 0.10, 
and ten cores have t/ tur buient/PI = 0.20 . 

3.2. Numerical method 

The simulations are performed with a standard SPH code 
(e.g. Lucy, 1977; Gingold & Monaghan, 1977; Monaghan, 
1992). For details of the code, see Turner et al. (1995) and 
Goodwin et al. (2004a,b). It uses variable smoothing lengths h, 
the M4 smoothing kernel (Monaghan & Lattanzio, 1985), and 
AA n eib =50 + 5 neighbours. The artificial viscosity prescription 
of Lattanzio et al. (1986) is used with a v i sc = 1 and /? v ; sc = 2. 
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Fig. 2. The stellar IMF produced from an ensemble of fifty core collapse simulations is shown as a solid-line histogram. 
This histogram is well-fitted by the sum (dashed line) of two distributions (dotted lines): a log-normal distribution with mean 
{ogio[M] = -0.05 and standard deviation o-( g m [M] = 0.02 plus a flat distribution below 0.5M Q The stars making up the log- 
normal distribution generally tend to remain in multiple systems, whilst those making up the flat distribution are generally objects 
that have been ejected from multiple systems. 



An octal tree (Barnes & Hut, 1986) is used to construct neigh- 
bour lists and to calculate gravitational accelerations, which are 
kernel softened using h. If the density of particle i rises above 
100p cl it, and the particles within 5 AU of particle i are gravita- 
tionally bound, they are replaced with a sink particle of radius 

5 AU (Bate, Bonnell & Price, 1995). 

Most simulations are made with N to t — 25,000 particles, 
but some have been repeated with 7V tot = 50, 000 particles in 
order to test convergence. In no case is there a significant differ- 
ence in the evolution using larger A/tot- With AAot = 25, 000, the 
minimum mass which can be resolved is M res 0.01 M Q . With 
the equation of state given by Eqn. (|2jl, the minimum Jeans 
mass is M m ; n =2 0.01M Q . Therefore the Jeans condition (Bate 

6 Burkert, 1997; Whitworth, 1998) is satisfied. 



3.3. Results 

We have previously shown (Goodwin et al., 2004a,b) that 
even cores with very low levels of turbulence can fragment to 
form multiple systems. In the ensemble of fifty simulations of 
Taurus-like cores reported here, each core forms between 1 and 
9 objects, with an average of 4.5 objects. Details of the individ- 
ual simulations are given in Goodwin et al. (2004a,b). Typically 
the lower-mass objects are ejected from the core in a period of 
rapid dynamical evolution (cf. Sterzik & Durisen, 2002), and 
the higher-mass objects remain in the core and form multiple 
systems. The ejected objects are lower-mass because they are 
ejected before they can grow much by accretion (cf. Reipurth & 
Clarke, 2001; Bate, Bonnell & Bromm, 2002; Delgado-Donate 
et al., 2003, 2004; Goodwin et al., 2004a,b). The objects re- 



maining in the core are higher-mass because they continue to 
accrete (cf. Bonnell et al., 2001). 

Figure|2]shows as a histogram the mass function of objects 
produced by this ensemble of runs. The histogram is well fitted 
by the sum (dashed line) of two distributions (dotted lines): 
a narrow log-normal distribution, with mean (iogio[M]) = 
-0.05 and standard deviation <T( ogl0 [M] - 0.02; and a flat distri- 
bution below 0.5M o . The stars contributing to the log-normal 
distribution are mainly those that remain in the core and end up 
in multiple systems. The stars contributing to the flat distribu- 
tion are mainly those ejected from the core. 

The low-mass objects have a flat mass distribution because 
their ejection is a stochastic process. In any individual core 
there are only a few low-mass objects, so the ejection dynam- 
ics cannot be very selective about which object is ejected (i.e. 
ejection is only weakly dependent on mass). The mass of an 
ejected object depends mainly on how long after formation it is 
ejected. If ejection occurs early, the object has had little time to 
grow and is likely to be a brown dwarf. If ejection occurs later, 
the object may have grown to become a low-mass star. 

4. Comparison of simulations with observations 

Figure [3] shows a comparison between the observed IMF of 
Taurus from Figure 1 (solid line histogram) and the model IMF 
from Figure 2, normalised and rebinned to the same binning as 
the observations (shaded histogram). The two histograms are 
very similar. A Kolmogorov-Smirnov test estimates the proba- 
bility that the underlying distributions are different at less than 
10%. Ay 2 test accepts the hypothesis that they are the same at 
more than 90% confidence. 
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Fig. 3. The solid-line histogram shows the observed IMF of Taurus due to Luhman et al. (2003a) (as in Figure 1). The shaded 
histogram shows the model IMF from Figure 2, normalised and rebinned to match the binning of Luhman et al. (2003a). The 
model IMF was produced from an ensemble of simulations which starts with the observed core mass distribution in Taurus and 
initial conditions for individual cores (density profile, level of turbulence, etc.) based on the physical conditions observed in 
pre-stellar cores in Taurus. The two histograms have been normalised by minimising y 1 . Agreement between the model and the 
observations is good (see text for discussion). 



The main discrepancy between the histograms is at high 
masses, where the simulations appear to produce too many 
stars above 1M . This may be due to the simulations not in- 
cluding feedback. Bipolar outflows from young stars would re- 
duce accretion onto the more massive stars, and might thereby 
bring the histograms into better agreement. We are currently 
repeating a subset of the simulations with a more sophisticated 
sink model which includes feedback due to bipolar outflows 
(Boyd et al., in preparation), in order to test this hypothesis. 

Briceno et al. (2002) report that the spatial distributions 
of stars and brown dwarfs in Taurus are very similar, and this 
would appear to be evidence against the hypothesis that brown 
dwarfs are formed by ejection from star forming cores (e.g. 
Reipurth & Clarke, 2001; see also Kroupa & Bouvier, 2003, 
for a discussion in the specific context of Taurus). A similar 
situation is seen in Chamaeleon I by Lopez Marti et al. (2003). 

However, in our simulations the ejected objects comprise 
both brown dwarfs and low-mass stars, with almost twice as 
many low-mass stars (0.08M o < M < 0.5M o ) as brown dwarfs 
(M < 0.08M G ). Moreover, the ejection speed is almost inde- 
pendent of mass (~ 1 - 2kms _1 ), and therefore the distribu- 
tions and kinematics of brown dwarfs and low-mass stars are 
very similar, as observed. 

Moreover, any systematic difference between the velocity 
dispersions of objects of different mass is unlikely to produce 
significant segregation on timescales less than 10 Myr, for the 
following reasons. First, the star-forming cores in Taurus are 
widely spaced, and so there is no single centre from which to 
measure a putative diaspora of brown dwarfs or low-mass stars. 
Second, the individual star-forming cores have different veloci- 



ties; from Onishi et al. (2002), the inter-core velocity dispersion 
is ~ ± 1 kms -1 . Third, some of the older, smaller star-forming 
cores may already have dispersed. 

5. Conclusions 

Taurus has an unusual stellar initial mass function and an un- 
usual core mass function. We have modelled the hydrodynami- 
cal evolution of an ensemble of cores with masses based on the 
Taurus core mass function and levels of turbulence based on 
those observed in Taurus. We find that the unusual stellar IMF 
in Taurus can be explained as a direct result of the unusual core 
mass function and intrinsic core properties in Taurus. 

In each core an initially unstable, multiple system forms, 
with between 2 and 9 members. Typically, 2 or 3 objects are 
ejected before they can accrete a significant amount of material. 
These ejected low-mass stars and brown dwarfs constitute the 
flat, low-mass 'tail' of the Taurus IMF. There are almost twice 
as many low-mass stars (0.08M© < M < Q.5M G ) as brown 
dwarfs (M < 0.08M G ). The ejection velocities (~ 1 - 2kms _1 ) 
are essentially independent of mass, so the low-mass stars are 
as dispersed as the brown dwarfs and there is no significant 
mass segregation. 

The remaining objects stay near the centre of the core and 
continue to accrete until their masses are ~ 0.8M o , by which 
stage there is not much mass left to accrete. These more mas- 
sive stars constitute the Gaussian peak centred at ~ 0.8M G in 
the Taurus IMF. Almost all of them are in binary or triple sys- 
tems. Thus the simulated cluster of cores has an IMF very sim- 
ilar to that of Taurus, viz. a narrow Gaussian peak at ~ 0.8M Q , 
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and a flat tail at lower masses, extending into the brown dwarf 
regime. 

Further support for the hypothesis that the form of the IMF 
in Taurus is a direct result of the typical masses of the cores 
in Taurus comes from Goodwin et al. (in preparation) who 
find that a core with low levels of turbulence typically forms a 
number of objects approximately equal to the number of initial 
Jeans masses in the core (the initial Jeans mass being roughly 
1M Q ). Thus lower-mass cores form fewer objects and are un- 
able to populate the tail of the IMF as they do not eject sig- 
nificant numbers of brown dwarfs and low-mass stars. Higher- 
mass cores over-populate the tail through the ejection of more 
brown dwarfs and low-mass stars. In addition the stars that re- 
main within a core have a larger resovoir of gas to accrete and 
become larger than the observed O.8M peak. 

The agreement between the two histograms of Figure 3 (ob- 
servation and theory) is remarkably close. Hence we believe we 
have shown, for the first time, a direct causal link between the 
core mass function of a star-forming region and the stellar IMF 
produced in that region. 
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